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Abstract. — Isoetes is an ancient lycopod lineage with a highly conserved morphology that 
provides few morphological characters to resolve the phylogeny of its species. Species appear to 
have evolved by divergence and allopolyploidy. The basic diploids I. hypsophilo, /. toiivanensis, 
and /. yunguiensis and the tetraploid /. sinensis occur in China. Analysis of ITS sequences 
indicates that the Chinese Isoetes species are part of an Australasian clade including /. brevicula 
from Western Australia and /. kirkii from New Zealand. Two distinct cloned sequences of the 
second intron of a homolog were recovered from /. sinensis supporting the hypothesis that 

/. sjnensis is an allotetraploid. One of the /. sinensis cloned sequences was similar to the /. 
taiwanensis sequence and the other cloned sequence was similar to the /. yunguiensis sequence 
identifying /. taiwanensis and I. yunguiensis as the likely parents of /. sinensis. Other cloned 
sequences recovered from I. sinensis were recombined parts of the two distinct sequences. 
Morphological evidence supporting an allotetraploid origin of /. sinensis is found in its larger 
microspore size and intermediate megaspore texture compared to 1. taiwanensis, and L 
yunguiensis. 


Isoetes L. is a cosmopolitan genus of heterosporous lycopods containing 
hundreds of species. Plants usually appear as tufts of linear leaves arising from 
an underground, corm-like rootstock. Ellipsoidal sporangia occur in expanded 
leaf bases. Species range from evergreen aquatics to ephemeral terrestrials. 
Although Isoetes is an ancient lineage with its distinctive morphology rec¬ 
ognizable in the Triassic (Retallack, 1997), few characters have been found in 
its highly conserved morphology to resolve the phylogenetic relationships of 
its species. Distinguished by their habitat preference, megaspore morphology, 
and chromosome numbers, species appear to have evolved by ecological 
isolation and genetic divergence as separated populations adapted to ter¬ 
restrial or aquatic habitats and by interspecific hybridization and chromosome 
doubling (allopolyploidy) when divergent species were dispersed into the 
same sites (Taylor et al., 1993). 

Interspecific hybridization and allopolyploidy are well documented for 
Isoetes. Many interspecific hybrids have been recognized by their production 
of irregular spores and confirmed by their chromosome numbers (Brunton and 
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Britton, 1999). In several cases, interspecific hybrids and their suspected 
allopolyploid derivatives have been verified by chromosome counts, isozyme 
profiles, and DNA sequences (Taylor and Hickey, 1992; Hoot and Taylor, 2001; 
Hoot et ah, 2004). A polyploid series ranging from 3.y = 33 to 12,y =132 is 
known for Isoetes. Over 60% of Isoetes taxa, for which chromosome counts 
have been published, are polyploid (Troia, 2001). Therefore, not only is there 
documentation that interspecific hybridization and allopolyploidy occur in 
Isoetes, but there is also evidence that allopolyploidy is an equally important 
mechanism of speciation in this genus. 

Recentlv, herbarium, field, and laboratory studies have been conducted to 
learn more about the Isoetes of China. These studies have provided an opportu¬ 
nity to determine the status of historical populations, discover new popida- 
tions and new taxa, and obtain live specimens from which root tips could be 
harvested for chromosome counts and fresh leaves could be collected for DNA 
isolations. 

Four species of Isoetes have been described for China. All are believed to be 
rare and endangered. Isoetes hypsophila Handel-Mazzetti is known from the 
Hengduan Mountains in northwestern Yunnan Province and southwestern 
Sichuan Province. In this region, /. hypsophila occurs in the shallow water of 
lakes and ponds about 3600 m above sea level. Isoetes sinensis T. C. Palmer has 
been found at about ten sites in and along rivers and lakes of the middle and 
lower Yangtze River system. At present, only three populations are known to 
remain in China. Isoetes sinensis has also been reported from the Kyushu and 
Chubu Districts Japan (Takamiya et al., 1997) and Cheju Island. South Korea 
(Takamiya, 2001). Isoetes taiwanensis DeVol is known only from Menghuan 
Lake near the foot of Zhixing Mountain in the Yangming Mountains National 
Park, north of Taipei in northern Taiwan. Isoetes yunguiensis Q. F. Wang and 
W. C. Taylor is known from the Yunnan-Guizhou Plateau in southwest China. 
In this region, plants have been recorded at four sites, but only two small 
populations, totaling about 400 individuals, are known to remain. Liu et al. 
(2002) reported that /. hypsophila, /. taiwanensis, and /. yunguiensis are basic 
diploids (2n = 2x = 22) and Isoetes sinensis is a tetraploid (2/7 = 4x = 44). 

Nucleotide sequences from the nuclear ribosomal ITS region, the chloroplast 
atpB-rbcL spacer region, and the second intron of a LEAFY homo\og have been 
used to determine phylogenic relationships of Isoetes, delimit species, and 
reveal an interspecific hybrid and its allotetraploid deriv^ative (Hoot and Taylor, 
2001). Hoot et al. (2004) used cloning to separate homoeologous sequences of 
the second intron of a LEAFY homolog for several Isoetes allotetraploids. By 
comparing these cloned sequences to those of putative parents, some of the 
parent species could be identified. 

The goals of this paper were to use nucleotide sequences from the 
nuclear ribosomal ITS region and the second intron of a LEAFY homolog 
to: (1) determine the relationships of the Chinese Isoetes species. (2) test the 
hypothesis that the tetraploid /. sinensis is an allotetraploid and, if this 
hypothesis is correct, (3) identify the basic diploid parent species of 
7. sinensis. 


198 


AMERICAN FERN JOURNAL: VOLUME 94 NUMBER 4 (2004) 


Table 1. Specimens sampled. Columns indicate species, location, collector, collection number- 
DNA isolation number, date of collection, and herbarium acronym for location of voucher. 
Collections are from Mainland China unless otherwise noted. 


Species 

Voucher Collection 

Isolation 

Number 

Isoetes brevicula 

Rock pool, summit of Lily McCarthy Rock, Western 

Australia, 25 Sep 2002, W. C. Taylor & N. T. Liiebke 

6383 (MIL) 

189 

Isoetes hypsophila 

Tu-er-sjan, Dao-cheng County, Sichuan Province, 01 Aug 
2001, Wang Aing-Feng, Liu Xing, Liu Hong & Yang 
Xiao^Lin 2 (WH) 

127 

Isoetes kirkii 

Lake Brunner, South Island, New Zealand, 27 Mar 2004, 

D. W. Woodland & Felicity Cutten s,n. (MIL) 

kiNZ 

Isoetes sinensis 

Xing-an-jiang, Jiande City, Zhejiang Province, 19 Oct 

2001, Liu Xing 8r Pang Xin-An 3, 4 (WH) 

129, 131 

Isoetes taiwanensis 

Menghuan Lake, Yangming Mountains, Taiwan, May 

1998, Cbiou Wen-Liang s.n. (MIL) 

78 

Isoetes yunguiensis 

Sha-shi-chong, Ping-ba County, Guizhou Province, 15 

Aug 2001, Liu Xing 8r Yang Xiao-Lin 5 (WH) 

130 


Materials and Methods 

Species sampling. —Table 1 contains locality, collector, collection number, 
and location of voucher specimens for plants used in this study. Specimens were 
identified to species using the original descriptions of the species (Handel- 
Mazzetti, 1923; Palmer, 1927; DeVol 1972a; Wang Q. F. et al, 2002) and by 
comparison with authentic and type specimens. Diagnostic morphology for 
megaspore textures was evaluated using an Olympus SZX12 stereomicroscope. 

DNA isolation and amplification. —DNA was isolated from 20 mg of silica 
dried leaves from each sample by grinding the leaf tissue, frozen in liquid 
nitrogen, to a powder with a disposable 1.5 pellet pestle (Kimble-Kontes) in 
a 1.5 ml snap-cap microcentrifuge tube (Eppendor^ and using the DNeasy*^ 
Plant Mini Kit (Qiagen) following the manufacturer’s protocol. 

The ITS region for all samples was amplified with the primers lTS-1 (5'- 

GTCCACTGAACCTTATCATTTAG-3'; Urbatsch, et al. 2000) and 1TS4 (5'- 
TGCTCCGGTTATTGATATGG-3'; White et al. 1990). The LEAFY intron for 
all samples was amplified with the primers 30F (5'- GATGTTTATGAA- 
CAATGTGG-3') and 1190R (5’- GAAATAGCTGATTTGTAACC-3'); Nancy S. 
Napier designed both LEAFY primers. PGR reaction mixtures followed 
manufacturer protocols using Ready-To-Go'*''' PGR Beads (Amersham Bio¬ 
sciences). PGR amplification for the ITS region began with denaturation for 60 
s at 97 C followed by 40 cycles of denaturation tor 10 s at 97'^G, annealing for 
30 s at 48 G, and extension for 20 s at 72°G with 4 s added to extension time 
each cycle and ending with a final extension of 7 min at 72°C. PGR 
amplification for the LEAFY intron began with denaturation for 5 min at 94°G 
followed by 40 cycles of denaturation for 1 min at 94°C, annealing for 1 min at 
50 C, and extension for 1 min at 72°C, and ending with a final extension of 5 
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min at 72°C. PCR products were concentrated via electrophoresis in a 2% 
agarose gel containing ethidium bromide and visualized with transilluminated 
UV. Bands were cut from the gel and purified using the QlAquick Gel 
Extraction Kit (Qiagen). 

Cloning and sequencing. —Sequencing of the ITS PCR products was 
performed on both 5' and 3' DNA strands using the amplification primers 
ITS-I and ITS4 as cited above. 

Ligation of the purified LEAFY PCR product and subsequent transformation, 
cloning, and visualization of transformed clones followed manufacturer 
protocols using the pGEM®-T Easy Vector Systems (Promega) with LB Amp 
100 X-gal plates (Teknova). 

For basic diploid species, /. tahvanensis and /. yunguiensis, eight clones 
(colonies) for each species were sequenced. For the tetraploid species, /. 
sinensis, 16 clones (eight from each of two plants) were sequenced to increase 
the odds of capturing all parental cloned sequences. For the outgroup species, 
/. brevicula E. R. L. Johnson, 1. hypsophila, and /. kirkii A, Braun, four clones 

from each species were sequenced. 

Cleaning and concentration of the vector DNA followed manufacturer 
protocols using the QIAPrep® Spin Miniprep Kit (Qiagen). Sequencing was 
performed on both 5' and 3' DNA strands using sequencing primers M13F and 
M13R in conjunction with the ABI Big Dye® Terminator Cycle v 3.1 Sequencing 
Kit (Applied Biosystems) following the manufacturer’s protocol. Sequencing 
products were resolved on an ABI (model 377) DNA sequencer at the Iowa State 
University DNA Sequencing and Synthesis Facility, Ames, Iowa. 

Sequence alignment and phylogenetic analysis. —Nucleotide sequences 
were aligned and edited using Sequencher 4.1 (Gene Codes Corp.). Gaps were 
treated as additional presence/absence characters, with one or multiple base 

gaps scored as a single character (Baldwin et al. 1995). 

Maximum parsimony analysis of the data was conducted with PAUP* 
version 4.0bl0 (Swofford, 2002) using the heuristic search option for the ITS 
data set and the LEAFY data set, maximum trees = 4000 for the ITS data set and 
maximum trees = 100 for the LEAFY data set. PAUP* was used to run 500 
bootstrap replicates tor each data set to estimate reliability of the clades 

(Felsenstein, 1995). 


Re.sults 

ITS sequences of Isoetes hypsophila, 1. tahvanensis and 1. yunguiensis from 
China. I. brevicula from southwestern Australia, and /. kirkii from New 
Zealand form an Australasian clade with other species and clades previously 
reported by Hoot and Taylor (2001). The data set analyzed consisted 16 
ingroup species and two outgroup species with a total of 826 characters; 313 
characters were variable and 239 were parsimony informative. The ITS tree 
illustrated is a bootstrap 50% majority-rule consensus tree of the 39 most 
parsimonious trees retained (Fig. 1). Tree topology shows a close relationship 
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Fig. 1. Isoetes ITS region tree. Bootstrap 50% majority-rule consensus tree of 39 trees resulting 
from maximum parsimony analysis using heuristic search of ITS region sequence data for eighteen 
basic diploid species of Isoetes. Tree length is 489 steps of equally weighted nucleotide 
substitutions and gaps, Cl (excluding uninformative characters) = 0.76, RI = 0.86. Numbers above 
the branches are the number of nucleotide substitutions. Numbers below the branches are 
bootstrap values. Major clades are labeled on the right. The tree is rooted with /. capensis and /. 
stellenbossiensis from South Africa. 


exists between /, taiwanensis, and I. yunguiensis. All species in the tree are 
basic diploids (2n = 22). 

Defined by the 30F and 1190R Napier primers, the LEAFY region for aligned 
sequence clones of I. sinensis, I. taiwanensis, and /. yunguiensis was 1105 
bases long and included all of the second intron and parts of the flanking 
exons. All of the eight sequence clones of /. taiwanensis were 1075 bases long 
and all of the eight sequence clones of /. yunguiensis were 1072 bases long. 
Two of the 16 /. sinensis clones did not amplify for the LEAFY intron. 
Therefore, 14 aligned, cloned sequences were compared for informative sites. 
At 15 informative sites that included 12 substitutions, two one base gaps, and 
one three base gap, six cloned sequences of /. sinensis, each 1076 bases long, 
matched or nearly matched the cloned sequence type of L taiwanensis and five 
cloned sequences of /. sinensis, each 1079 bases long, matched or nearly 
matched the consensus cloned sequence type of L yunguiensis (Table 2). Seven 
/. sinensis cloned sequences showed evidence of recombination i.e., part of the 
PCR amplified, cloned sequence first matched either the /. taiwanensis or the 
L yunguiensis sequence type, but further on matched the other sequence type. 
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Table 2. Comparison of fourteen cloned sequences of /. sinensis with cloned sequences of 
I. taiivanensis and I. yunguiensis at fifteen informative sites. Sites are numbered sequentially 
beginning with the 30F Napier primer. Isoetes sinensis cloned sequences are identified by 
collection number-DNA isolation number-clone number. Nucleotides in italic compare to those of 
/. taiwanensis. Nucleotides in bold face compare to those of /. yunguiensis. Gaps are indicated by 
dashes. Sequences marked with an asterisk showed extensive recombination and were removed 
from the data set before the final analysis. Clones 4-131-6 and 4-131-7 did not amplify for the 
LEAFY region and are not included in the table. 


site/clone 

205 

289 

290 

357 

525 

549 

660 

741 

776 

827 

840 

865 

965-967 

987 

1059 

taiwanensis 

c 

G 

T 

T 

G 

G 

T 

T 

A 

A 

A 

A 

• - - 

- 

A 

3-129-3 

c 

G 

T 

T 

G 

G 

T 

T 

A 

A 

A 

A 

• - - 

• 

A 

3-129-7 

c 

G 

T 

T 

G 

G 

T 

T 


A 

A 

A 

. • • 

• 

A 

4-131-1 

c 

G 

T 

T 

G 

G 

T 

T 

A 

A 

A 

A 

• • - 

• 

A 

4-131-3 

c 

G 

T 

T 

G 

G 

T 

T 

A 

A 

A 

A 

... 

• 

A 

4-131-4 

c 

G 

T 

T 

G 

G 

T 

T 

A 

A 

A 

A 

• • • 

• 

A 

3-129-5 

c 

G 

T 

T 

C 

G 

T 

C 

A 

A 

A 

A 

• - - 

• 

G 

3-129-2* 

c 

G 

T 

T 

T 

T 

C 

c 

G 

G 

- 

G 

• • . 

- 

A 

3-129-1* 

c 

G 

T 

T 

T 

T 

C 

c 

G 

G 

- 

G 

A T G 

T 

A 

4-131-8* 

T 

A 

C 

C 

C 

G 

T 

T 

A 

A 

A 

A 

- - - 

• 


4-131-2 

c 

A 

C 

c 

T 

T 

C 

c 

G 

G 

- 

G 

A T G 

- 

G 

3-129-4 

T 

A 

c 

c 

T 

T 

C 

c 

G 

G 

- 

G 

A T G 

- 

G 

3-129-8 

T 

A 

c 

c 

T 

T 

c 

c 

G 

G 

- 

G 

A T G 

• 

G 

3-129-6 

T 

A 

c 

c 

T 

T 

c 

c 

G 

G 

- 

G 

A T G 

T 

G 

4-131-5 

T 

A 

c 

c 

T 

T 

c 

c 

G 

G 

- 

G 

A T G 

T 

G 

yninguiensis 

T 

A 

c 

c 

T 

T 

c 

c 

G 

G 

- 

G 

A T G 

T 

G 


Three sequences showing evidence of extensive recombination were removed 
from the data set evaluated by PAUP*. 

Based on the relationships indicated in the ITS tree (Fig. 1), /. brevicula, I. 
kirkii, and I. hypsophila were chosen as outgroup species to root the LEAFY 
second intron tree. The sequenced LEAFY region for /. brevicula was 1035 
bases, for /. kirkii it was 1090 bases, and for /. hypsophila it was 1095 bases. 

The LEAFY second intron data set analyzed contained 16 sequences with 
a total of 1125 characters: 171 characters were variable and 54 were parsimony 
informative. The data set included consensus sequences from clones of /. 
brevicula, I. hypsophila, 1 . kirkii, /. taiwanensis, and I. yunguiensis and 11 cloned 
sequences of /. sinensis. The LEAFTsecond intron tree illustrated is a bootstrap 
50% majority-rule consensus tree of the two most parsimonious trees retained 
(Fig. 2). Six of the/, s/nensis cloned sequences formed a clade with /. taiwanensis 
and five of the /. sinensis cloned sequences formed a clade with /. yunguiensis. 

Discussion 

Although Isoetes is Paleozoic in origin (Pigg. 1992), worldwide in 
distribution, and over time, undoubtedly adapted to changing climates and 
aquatic to terrestrial habitats on every continent many times, the morphology 
of Isoetes has been remarkably conserved. Thus, morphology provides few 
characters that can be used to reliably reconstruct phylogenetic relationships. 
Nevertheless, pteridologists have speculated about the relationships of Isoetes 
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Fig. 2. Isoetes LEAFY' second intron Iiomolog tree. Bootstrap 50% majority-rule consensus tree of 
four trees resulting from maximum parsimony analysis using lieuristic search of LEAFY' second 
intron homolog sequence data for five basic diploid and one tetraploid [I. sinensis) species of 
Isoetes. Eleven sequence clones that form a sister clade to either /. yunguiensis or 1 . taivvanensis 
represent the tetraploid genome of /. sinensis. Tree length is 183 steps of equally weighted 
nucleotide substitutions and gaps. Cl (excluding uninformative characters) = 0.89, RI = 0.96. 
Numbers above the branches are the number of nucleotide substitutions. Numbers below the 
branches are bootstrap values. Figures to the right of the specific epithets are the collection and 
clone identification labels. The tree is rooted with /. hypsophila from China, /. breviculo from 
Western Australia, and /. kirkii from New Zealand. 

species based on ecology, morphology, and biogeography. Britton and Brunton 
(1991) reevaluated the spore morphology of /. tahvanensis, concluding that it 
was not related to taxa from southwestern Australia as proposed by Marsden 
(1979), but instead appeared to have its closest affinity to /. kirkii from New 
Zealand. The ITS tree (Fig. 1) shows that both /. breviculo from southwestern 
Australia and /. kirkii form a sister clade to the Chinese species and all are 
members of an Australasian clade. Based on spore morphology, habit, and 
habitat, Huang et al. (1992) concluded that /. taiwanensis is probably closer to 
/. asiatica than it is to /. sinensis, but /. asiatica (/. echinospora subsp. asiatica 
(Makino) A. Love is a member of the I. echinospora species complex, a group of 
circumpolar taxa with echinate megaspores (Love, 1962, Takamiya, 1997). The 
ITS tree (Fig. 2) shows that /. echinospora, a member of an American clade, is 
only distantly related to members of the Australasian clade, which includes /. 
taiwanensis. DeVol (1972h) mentioned that /. taiwanensis seemed nearer to 1. 
sinensis than any other species. Takamiya (2001) saw similarities in the spore 

morphology of /. taiwanensis and /. sinensis and concluded that phylogenetic 
comparisons of these two taxa were needed. 
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Fig. 3. Hypothetical phylogeny of Isoetes sinensis involving interspecific hybridization and 
chromosome doubling of the basic diploids /. taiwanensis and I. yunguiensis. Data presented in 
this study supports the hypothesis of an allotetraploid origin for /. sinensis. 


The recovery of two, distinct, LEAFY second intron sequence types from the 
tetraploid /. sinensis supports the hypothesis that /. sinensis is an allotetraploid 
(Table 2). To clarify the results, it was assumed that some of the sequenced 
clones recovered from /. sinensis were recombinations of the two distinct 
sequence types and three of the recombinant sequences were removed from the 
data set evaluated by PAUP*. Since recombination occurs from crossing over 
between chromosomes during meiosis, it is possible that the observed re¬ 
combined sequences were products of natural events. If the recombined 
sequences detected were the result of crossing-over during meiosis, we would 
predict that identical crossover sequences would be recovered as clones. All 
seven of the recombinant cloned sequences from I. sinensis were different, 
indicating that these recombined sequences more likely occurred during the 
PCR amplification reaction. Whatever their source, including recombinant 
sequences in a cladistic analvsis will affect results and therefore, thev need to be 
recognized and removed from the data set before the final analysis. 

Comparison of the two distinct cloned sequences from the allotetraploid /. 
sinensis with the cloned sequences of /. taiwanensis and /. yunguiensis indi¬ 
cates that either these two basic diploid species, or closely related taxa, likely 
participated in the formation of /. sinensis (Table 2). Although the /. 
yunguiensis clade, including five /. sinensis cloned sequences, and the /. 
taiwanensis clade, including six I. sinensis cloned sequences, are both well 
supported with high bootstrap percentages, the /. yunguiensis sequence is 
distinguished from its sister /. sinensis clones by sixteen autapomorphies and 
the /. taiivanensis sequence is distinguished from its sister I. sinensis clones by 
two autapomorphies (Fig. 2). These unique nucleotide substitutions could be 
due to (1) sequencing nucleotides of ta.xa different from those of the parent taxa, 
(2) continued evolution of the progenitor parent species and the allotetraploid 
species following allopolyploidy, or (3) copy errors during PCR. Causes for the 
autapomorphies might be determined by additional sampling and repeated 
PCR of the same clones. 

In addition to the molecular characters, morphological characters indicate 
that /. sinensis, I. taiwanensis, and /. vuneuiensis are distinct, but closely 
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related species and provide some evidence supporting the allopolyploid origin 
of /. sinensis. All three species are amphibious plants with tri-lobed rootstocks. 
They all have a rudimentary velum that covers only the upper edge of the 
sporangium. The microspores of /. taiwanensis and /. yunguiensis range from 
20-26 pm in length whereas, those of /. sinensis range from 26-30 pm in length 
(Britton and Brunton, 1991; Wang et al., 2002; Palmer, 1927). The larger size of 
/. sinensis microspores is attributed to its increased chromosome number. 
Increases in chromosome number are usually accompanied by larger spore size 
(Kott and Britton, 1983). In contrast, megaspores of the tetraploid /. sinensis 
and the basic diploid I. yunguiensis average about 400 pm in diameter, whereas 
megaspores of the basic diploid /. taiivanensis average about 300 pm in 
diameter. Megaspore texture appears to be the most distinctive character that 
separates these three species. However, in view of the results presented here, 
the cristate to verrucate megaspores of I. sinensis can also be interpreted as 
subtly combining the textures of the rugulate to reticulate megaspores of /. 
taiwanensis and the cristate to reticulate megaspores of /. yunguiensis. 
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Appendix 

Genbank accession numbers for the new Isoetes DNA sequences used in this manuscript. 


Species 

Genbank 

accession number 

Nuclear ribosomal ITS region sequences 

/. brevicula 

AY641098 

I. hypsophilo 

AY641099 

I. kirkii 

AY641100 

/. taiwanensis 

AY641101 

I. yunguiensis 

AY641102 

LEAFY second intron homolog sequences 

I. brevicula 

AY64tl03 

I. hypsophilo 

AY641104 

I. kirkii 

AY641105 

I. taivvonensis 

AY641106 

I. yunguiensis 

AY641107 

I. sinensis (/. yunguiensis type clone) 

AY641108 

I. sinensis (I. taivvanensis type clone) 

AY641109 







